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The work extractable from correlated bipartite quantum systems can be used to distinguish en- 
tanglement from classical correlation. A natural question is now whether it can be generalized to 
multipartite systems. In this paper, we devise a protocol to distinguish the GHZ, the W, and sepa- 
rable states in terms of the thermodynamically extractable work under local operations and classical 
■ communication, and compare the results with those obtained from Mermin's inequalities. 

<N ; 

! I. INTRODUCTION 

OH' 

Finding an efficient separability criterion for multipartite quantum systems and characterizing entanglement have 
been an important problem in the field of quantum information theory. Since Bell formulated the statistical irregular- 
ities that can be seen in the correlations between measurement outcomes from two distant, but entangled, systems 
various criteria and related notions have been discovered in the context of non-local correlations in quantum systems 
0, S 0, HI • Recently, it has been shown that it is possible to devise a protocol to extract more work via entangled 
qq ' systems from a heat bath thermodynamically than can be done from any separable state [||. This is of interest 
because a physically useful quantity, locally extractable work, can be employed to test the existence of entanglement 
t-H ' in bipartite quantum systems, leading to the idea of thermodynamical separability criteria. 

In this paper, we discuss if the thermodynamical separability criteria can be generalized to multipartite systems, 
particularly tripartite ones, as even the simplest transition from bi- to tripartite systems makes our problems much 
harder. It has been known that, unlike bipartite systems, there are two non-equivalent classes of entanglement, i.e. 
the GHZ and the W states, when three quantum subsystems have non-local correlations 0. Suppose that we are 
given an ensemble of tripartite systems, which is in either the separable or the GHZ or the W states, where the GHZ 
state || can be in general written as 



in 
o 



|GHZ} = -±=(|000) + |111)), (1) 
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and the W state 0, 



|W) = -^(|001) + |010) + |100)), (2) 



where |0) and 1 1) are two orthogonal states of a two-level system, e.g. eigenstates of a Pauli spin operator a z . Our task 
here is to distinguish each of these two ensembles from separable states in terms of extractable work. By extractable 
work, we mean work that can be extracted thermodynamically from local heat baths under local operations and 
classical communication (LOCC). 



II. THE PROTOCOL FOR BIPARTITE CORRELATIONS 



Before discussing the case of tripartite entanglement, we here sketch the protocol shown in [6j for the work-extraction 
from bipartite quantum systems. If we have an ensemble of two-level systems, either classical or quantum, we can 
extract work of amount of fcTln2[l — H(X)\ from a heat reservoir thermodynamically |gL IToL ITU, where k,T, and 
H(X) are the Boltzmann constant, the temperature of the heat bath, and H(X) is the Shannon entropy of a binary 
random variable X, respectively. The variable X corresponds to the outcome of measurement on the system and 
H(X) can be written as H(X) — — po log 2 po ~ Pi log 2 Pi, where po and p\ are the probabilities for the two outcomes. 
We will set kT In 2 = 1 for simplicity hereafter and call the unit of work "bit" . We will also write " extracting work 
from quantum systems" for short, instead of "extracting work thermodynamically from a heat bath via quantum 
systems" . 
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FIG. 1: Schematic view of the protocol to extract work from correlated pairs. Dividing the whole ensemble into groups of 
two pairs, Alice and Bob use A(9) and B(6') for their measurement and work extraction. For a half of this ensemble, Alice 
measures her state with A(9) and Bob maximizes the extractable work from his side along the direction of 6' by using Alice's 
measurement results. For the other half, they exchange their roles. 



The extractable work from correlated pairs in [y] is a simple generalization of the above case. Suppose that two 
distant parties, Alice and Bob, have an ensemble of identically prepared pairs of quantum bits {qubits), which is 
described by a density operator p. First, we define the extractable work £ p (A(0), B(8')) when Alice and Bob chose 9 
and 9' as the directions of their (projective) measurement. After dividing the shared ensemble into groups of two pairs, 
Alice measures one of the two qubits in a group with the projector she chose and informs Bob of the outcome (See 
Fig. |2J). Bob performs the same on his qubit of the other pair in the group. As a result of collective manipulations, 
they can extract £ p (A(O),B(0')) = 1 — 1/2 • [H{A(9)\B{9')) - H{B{9')\A(9))] bits of work per pair at maximum. 

Second, we consider a quantity S(p), which is an average work extractable when we set 9 = 9' and vary 9 continuously 
over a great circle on the Bloch sphere: 

S(p):=l-J\ p (A(9) ) B(9))d9, (3) 

The great circle is the one that maximizes the integral. The integral in Eq. J3J can be taken over the whole Bloch 
sphere, nevertheless, it does not change the essential part of our discussion here. 
In ly], it was shown that an inequality 

E(p) < E(|00)) = 0.4427 (4) 

is a necessary condition for a two-dimensional bipartite state p to be separable, that is of the form, p = J2 t Pipf <8> pf ■ 
The state |00) in the right-hand side (RHS) can be any pure product state \ipip')- 

If we integrate £ p over the whole surface of the Bloch sphere, then the condition 0} becomes 

Zbs(p) < S(|00» = 0.2787, (5) 

where the subscript BS stands for the Bloch sphere. This condition, Eq. (J3J), has been shown to be more effective 
than the standard Bcll-Clauser-Horne-Shimony-Holt (Bell-CHSH) inequalities 0] m detecting the inseparability of 
the Werner state which is a state in the form of pw = p\^f~)(^f~\ + (1 — p)/4 ■ I- 



III. PROTOCOL FOR TRIPARTITE SYSTEMS 



Let us now consider the case of tripartite quantum systems. By generalizing the protocol in the previous section 
in a simple manner, we can have a necessary condition for separability, which is of the same form as Eq. 0}, 
straightforwardly. That is, each one of three parties receives the outcomes of measurements from the other two and 
we take conditional entropies such as H(A(9)\B(9), C(9)) instead of H(A(9)\B(6)) in the definition of £ and S. Then, 
the inequality 0} with S(|000)) on the left-hand side (and a different numerical value)holds for all separable states 
of the form of p = J2iPiPt ® pf ® Pi ' ■ The proof is essentially the same as that in [||. Thus any violation of the 
inequality implies the existence of entanglement between at least two subsystems. Obviously, this extension can be 
applied to correlations between larger number of subsystems. 

However, such an inequality tells little about the properties of multipartite entanglement. Thus, we would like to 
find a different protocol and focus on distinguishing the two inequivalent classes of tripartite entanglement: in this 
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FIG. 2: Work extraction from a tripartite system. Alice and Bob measure their particle along the direction of z and u, 
respectively, where Bob varies the direction of u, choosing one of the three orthogonal directions {x, y, z} at each round. 
Charlie extracts work along u after receiving measurement results from Alice and Bob. They exchange their roles cyclically 
after three rounds, however, we do not consider the exchange of roles as the states of our interest here are symmetric with 
respect to three parties. 



paper, we will devise a way to distinguish the GHZ state, the W state, and separable states, in terms of the extractable 
work. 

With three parties (Alice, Bob, and Charlie), we consider a protocol in which one extracts work along the direction 
of the it-axis after receiving information on the outcomes of the other two's measurements along the z- and u- axes, 
respectively. Namely, for a subensemble where Alice and Bob measure their qubits and Charlie extracts work from 
his, the (average) amount of work obtainable is 

wg )it (j>) = 1 - H(C(u)\A(z),B(u)). (6) 

In this equation, we denote the direction of axis by a unit vector, such as z. We have chosen two directions, z and it, 
for two measurements in anticipation that the difference between the GHZ and the W states will be seen by varying 
the direction of u, while keeping z fixed. We do not lose generality by choosing the direction of the work-extraction 
to be the same as one of the measurement axes, that is u in Eq. Jo}, as it is the right choice in detecting singlet-type 
bipartite entanglement. Also, as we focus on the GHZ and the W states, which are symmetric in terms of A, B, and 
C, we will assume throughout the paper that it is Charlie who extracts work after Alice and Bob make measurements 
on their subsystems. Otherwise, we need to permute the roles of each party and take an average. 

We now show that a quantity below, W(p), achieves our task. We assume that the direction z is given and fixed 
for everyone of three parties throughout our discussion. Instead of varying u continuously, we average the extractable 
work over three orthogonal directions in the Bloch sphere, (x,y,z), for u. This is mainly because of our limited 
resource for numerical computation, but it does not change our main task at all. With only the direction of z fixed, 
those of x and y cannot be determined. We use (j) to specify the angle between x and a certain predetermined direction 
in the space. Then, the averaged work can be written as 

and we take the maximum value of W^(p) over (j) to remove the ^-dependence as 

W(p)=maxW lj> (p). (8) 

We will analyse how we can use W(p) to distinguish three classes of correlations in tripartite systems. 

A. Product states 

With a product state, \^p A ^p B tp c }, the extractable work depends only on the subsystem from which we extract work, 
regardless of the measurements on the other two. Thus, W(\ip A ip B ip c }) will be determined solely by the geometric 
relation between the state of the system and the choice of axes. W c i ) (\ip A ip B ip c )) takes its maximum value of 1/3 
when the Bloch vector representing the state coincides with one of the axes, i.e. W(\ip A ip B ip c )) < 1/3 for a given z. 
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B. Separable states 



We now show that 1/3 is indeed an upper bound for the extractable work W from any separable state. That is, an 
inequality, 



W( P ) < \, 



(9) 



is a necessary condition for a tripartite state p to be separable. This condition is very similar to the thermodynamical 
separability criterion for bipartite correlations 0, however, W(p) is an average of work over only three directions 
as in Eq. J7J, while it was over a great circle on the Bloch sphere (or the whole surface of the sphere) in This 
results in a difference in the efficiency of detecting weakly entangled states: The condition JjJ) is less efficient than 
the criterion in [6j, in this sense. 

Proof. The conditional entropy in Eq. © can be written as 



H(C(u)\A(z),B(u)) = J2 K4» 4)^(^(^)14. 4) 

i,j = {0,l} 



(10) 



where At (Bi) means that Alice (Bob) obtained the outcome i (j) by the measurement along z (u) and p{At, Bi) is 
the probability of a joint event of At and Bi. Omitting the directions of axes for simplicity, we can write Eq. (|10|) as 



H(C\A,B) = Y,p ij H(C\A i ,Bj). 



fill 



If the state p is separable, i.e. p — ^ n PnPn ® Pn ® Pni then the density operator for Charlie after Alice and Bob 
obtained i and j becomes 

Pi = -T EP« Tr K^) Tr KPn )Pn, (12) 

Pi 7 
J n 

where PI {Pi) is a projection operator for outcome i (j) along the direction of z (u) and pij — 
EnPrMPzPn)Tr(PiPn)- Therefore, 

H(C\A,B) = J2 H (— E^ Tr ( P >«) Tr ( P >«) Tr ( P >»)) 

> J2IlP^( P ^nmPiPn)H(Pn) 

n i,j 

= Y^PnH(Pn), (13) 

n 

because of the concavity of the Shannon entropy. If we were sure that we had a tripartite product state, p^ ® p^ <£> p% , 
the conditional entropy would be 



1 



H(C«\A«,B«) = H ^Tr(PX)Tr(P^f )p° 



(14) 



where p*. = Tr(P>£)Tr(P^). 

Since W{p k ) = max[l - (1/3) J2a H(C k \A k , B k )} < 1/3 for all separable states p k , combining Eqs. ^ and fity. 
we have 



W(p) = max 



< max 



1 

* 3' 



l-\Y,H(C\A,B) 

U 

1 -l J E J EPnH(C»\A»,B») 



(15) 



thus Eq. @. As W(p) still has a dependence on the direction of z, even a highly entangled state may not violate the 
inequality 10, however, once it is violated it is surely a manifestation of entanglement. 
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C. GHZ and W states 

Let us discuss how we can distinguish the GHZ, and the W states, using the above criterion. To this end, we make 
use of the remaining variable, i.e. the direction of z. 

The extractable work from a GHZ state can reach 1 bit, which is the maximum possible violation of the inequality 
©. This maximum is attained when we choose z to be perpendicular to |0) in the Bloch sphere (such as |+) = 
(|0) + \l))/\^2), the state after Alice's measurement still has maximum (bipartite) entanglement. Because of the 
perfect correlation between B and C, Charlie can extract maximum work, 1 bit per one set of three subsystems. 

We have obtained the value of the minimum extractable work from the GHZ as 0.1619 bits numerically. This occurs 
when 9 = 1.5560 and <p = 0.5600, for example, where 9 and 4> are the azimuthal and the longitudinal angles in the 
Bloch sphere for the direction of Alice's measurement basis z. That is, the direction denoted by 9 = is that of |0) 
and <f> is a relative phase between |0) and |1) in their superpositions. 

The maximum work extractable from |W) is 7/9, which is attained when the indirection is parallel to |0). If the 
outcome of Alice's measurement is 0, which occurs with probability 2/3, the other two systems, B and C, will still 
have a perfect correlation along the three orthogonal directions, being in the state (l/v / 2)(|01) + 1 10) ) and thus Charlie 
can get 1 bit of work. If Alice's outcome is 1, then Charlie can obtain only 1/3 bits on average as the remaining two 
systems are in a product state |11). Therefore, the average extractable work will be (2/3) • 1 + (1/3) • (1/3) = 7/9 bits. 
The minimum work from the W state was obtained numerically as 0.1696 bits. This is achieved when 9 = 0.7169 and 
(j> = 7r/4 for the direction of z. 

Now we can distinguish the GHZ and the W states, looking at the maximum and the minimum values of W after 
varying the direction of z. The maximum and minimum extractable work from the GHZ and the W states have the 
following relationships: 

I < maxPF(|W)) < maxiy(|GHZ)), 

Q Z Z 

and minVK(IGHZ)) < minW(|W)), (16) 

z z 

although the difference between the two minima is rather small. Therefore, when the state of a given ensemble is one 
of the three possibilities, the GHZ, the W, and separable states, as we have assumed above, it is possible to specify 
the state by examining the range of extractable work. If W(p) < 1/3 always holds regardless of the direction of z, 
then it is in a separable state due to the above proposition. 

Suppose that entanglement exists only in two subsystems out of three, such as a = a AB ®a c , where a AB = |5'~)(\E' _ | 
is a maximally entangled state. If we wish to discriminate such bipartite entanglement in tripartite systems, we need 
to examine the extractable work from each party, instead of simply computing the average W. In this example of 
state (j, Alice and Bob can obtain the maximum work (1 bit), regardless of the outcome of Charlie's measurement, 
because of the perfect correlation between A and B. However, no information from Alice and Bob can be useful 
to maximize the work Charlie can extract. Thus, he can have only 1/3 bits of work at maximum after averaging 
over three directions. If one of the subsystems is disentangled from the other two, the extractable work from this site 
cannot exceed 1/3, which is precisely what has been proved in the preceding subsection. By comparing the extractable 
work from each site, we can distinguish tripartite correlation with bipartite entanglement. 

As we have mentioned above, the thermodynamical separability criterion for bipartite systems is able to detect 
more inseparability of the Werner state of two qubits than the Bell-CHSH inequalities. Let us consider a Werner-type 
state, in analogy of the bipartite Werner state, with tripartite systems [l4|: We define the GHZ- Werner- type state as 

/9 G HZ =p|GHZ)(GHZ| + i^./, (17) 

and the W- Werner-type state as 

^=p|W)(W| + ^-/. (18) 

These states are also referred to as isotropic states [l^. The inequality © turns out to be violated when p > 0.6521 
with the GHZ- Werner-type state, and when p > 0.6981 with the W- Werner-type state. Are these criteria better than 
other schemes of multipartite-entanglement-detection? In [l4(, the GHZ- Werner-type state with p > 0.3226 has been 
shown to be distillable, thus entangled, where p = (8/ - l)/7 with / = (GHZ|p^ HZ |GHZ) in [lj, while we have no 
corresponding data for the W- Werner-type state. 

Let us take Mermin's inequality as the multipartite version of Bell-CHSH inequalities Mermin's inequality for 
a tripartite system can be written as 



(B 3 ) = Tr(pB 3 ) < 2, 



(19) 
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TABLE I: The comparison of the thermodynamical separability criterion for tripartite systems, Eq. and Mermin's inequal- 
ity, Eq. dU, in detecting the inseparability of the Werner-type states. This table shows the minimum values of p to violate 
the inequalities, i.e. the smaller value, the more detection of inseparability. 





Therm, sep. criterion J3J 


Mermin's ineq. 119H 


GHZ- Werner-type state 
W- Werner-type state 


p > 0.6521 
p > 0.6981 


p > 0.5 
p > 0.6566 



where the operator B3 is defined as (by omitting the ® sign) 

"3 = (Oi CT 2 ' + a V a 2 )°3 + ( 1 2 - 1' 2' ) 3' ' ( 20 ) 

In Eq. I|20[) . af = at ■ a x is a measurement operator at site x € {A, B,C} and is a unit vector representing the 
direction of the measurement. 

The GHZ state, Eq. (Q, can violate the inequality l|19|) maximally as (B^)qhz = 4. This means that the GHZ- 
Werner-type state, Eq. 1171) . can violate Eq. I|19[) when p > 1/2 since 

(£3) GHZ- Weruer = Tr(p™ Z 6 3 ) = p(S 3 >GHZ- (21) 

Similarly, the W- Werner-type state, Eq. (|18ll . can violate the inequality when p > 0.6566 as the maximum value of 
(^3)w-Wcrner is 3.046 ^3] ■ Therefore, the condition © is less effective in detecting the entanglement in both the 
GHZ- Werner-type and the W- Werner-type states: There are a class of weakly entangled states that violate Mermin's 
inequality, but do not violate Eq. ©. 

The situation does not change much even if we average the extractablc work by varying the u over the whole Bloch 
sphere. The inequality © can be violated by the GHZ- Werner-type state when p > 0.8392 and by the W- Werner-type 
state when p > 0.9057, after choosing z optimally. The reason why Eq. © is not as effective as Mermin's inequality 
is yet unclear, while a similar thermodynamical criterion is more effective in the case of bipartite entanglement as 
shown in 0. 

IV. SUMMARY 

We have generalized the thermodynamical separability criterion to tripartite quantum systems and shown that 
it can be used to distinguish two different classes of tripartite entanglement, the GHZ and the W, in terms of 
thermodynamically extractable work under LOCC. We have also found that the criterion for tripartite systems is 
less effective in detecting the Werner- type entanglement than that for bipartite systems. Although it is not perfectly 
clear if the Werner-type states above are the proper counterpart of the bipartite Werner state we should compare 
with, it appears that the separability criterion we have obtained here is not as efficient as other criteria, such as 
Mermin's inequality, in detecting weak entanglement. Nevertheless, it is capable of distinguishing the GHZ, the W, 
and separable states with physical quantity and this may lead to an interesting physical or information processing 
process, in which the W state is more useful than the GHZ state, unlike most of the known processes. 
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